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g ' Abstract 

' We investigate the phenomenon of the decay of a supercurrent through ho- 

mogeneous nucleation of vortex-antivortex pairs in a 2-D like superconductor 
or superfluid by means of a quantum electrodynamic formulation for the de- 
cay of the 2-D vacuum. The case in which both externally-driven current 
^ | and Magnus force are present is treated exactly, taking the vortex activation 

£^ ' energy and its inertial mass as independent parameters. Quantum dissipation 

is included through the formulation introduced by Caldeira and Leggett. The 
most relevant consequence of quantum dissipation is the elimination of the 



threshold for vortex production due to the Magnus force. In the dissipation- 
dominated case, corresponding formally to the limit of zero inertial mass, an 
exact formula for the pair production rate is given. If however the inertial 
mass is strictly zero we find that vortex production is inhibited by a quan- 
tum effect related to the Magnus force. The possibility of including vortex 
pinning is investigated by means of an effective harmonic potential. While an 
additional term in the vortex activation energy can account for the effect of a 
finite barrier in the direction perpendicular to the current, pinning along the 
current depresses the role of the Magnus force in the dissipation-dominated 
dynamics, except for the above-mentioned quantum effect. A possible de- 
scription of vortex nucleation due to the combined effects of temperature and 
externally-driven currents is also presented along with an evaluation of the 
resulting voltage drop. 
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I. INTRODUCTION 



The layered structure and high critical temperature of the oxide supercon- 
ductors has lead to much renewed interest in the behaviour of these materials 
under applied magnetic fields. Of particular interest are the statics and dy- 
namics of the magnetic vortices in the presence of pinning impurity centers, 
where issues such as the formation of a glassy vortex phase M and those 
connected to transport properties such as flux-creep Q and vortex tunneling 
||-||] have been vivaciously debated in the recent literature. Similar questions 
of vortex dynamics have attracted the interest of researchers in the field of 
superfluidity for the past three decades |J|7)- 

An issue that appears to lend itself to full analytic treatment by means 
of a path-integral formulation is the quantum tunneling of a vortex trapped 
by a pinning center in a 2-D superfluid or superconductor structure 0,0]. 
This problem has been investigated both with and without the inclusion of 
quantum dissipation and discussions that give full weight to the inertial mass 
of the moving vortex can be found in the literature alongside treatments 
H in which the inertial mass is taken to be negligible. The physics of the 
problem is captured by semi-classical evaluations of the tunneling rate and 
discussion has centered on the effects of dissipation and pinning in opposing 
the stabilising influence of the Magnus force on the classical orbits of a vortex. 
Mathematically, and in the classical framework for a single vortex moving in 
a supercurrent at relatively low velocity, the equation of motion in the 2-D 
plane reads: 

mq = — VT/(q) — eqxB — ??q (1-1) 

Here m is the inertial mass of the vortex of topological charge e = ±2ir, 
treated as an single, point-like particle of 2-D coordinate q(i). Also, V(q) = 
^p(q)+K>(q> J) is the phenomenological pinning potential plus the electric-like 
potential due the supercurrent J, namely — W v = eE = xeJ (the latter being 
a compact notation for the 2-D dual of a vector: (xj)^ = e^Jj,). Finally 
B = zdp^ is the magnetic-like field of the Magnus force [p|,|8|-|lcfl acting on 



the vortex in the plane of the film of thickness d and directed along the 
vector z orthogonal to it and associated with a superfluid component having 

(3) 

3-D number density p s , and r/ is the phenomenological friction coefficient. 
We neglect in this work the additional magnetic field intrinsic to the vortex 



since, as discussed by Clem |11], the flux quantum is homogeneously spread 
over a very large distance for a single flux vortex in a 2-D superconductor. 
The quantum-mechanical version of Eq. ([O]) is attained via the Feynman 
path-integral transposition in which the dissipation is treated quantistically 
through the formulation due (for the ohmic case) to Caldeira and Leggett 



\m- 



A situation closely related to the quantum tunneling of a vortex in the 
presence of dissipation is the decay of a supercurrent, whether in a charged su- 
perconductor or in a neutral superfluid, due to the spontaneous homogeneous 
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nucleation of vortex-antivortex pairs from the "vacuum" of the perfect super- 
fluid wavefunction, in the presence of an external uniform field. In this case 
the field corresponds - in its electric part - to the dual of the external super- 
current J, which may be kept fixed in the system if a voltage drop is measured 
at the edge of the superconducting sample, and - in its magnetic part - to the 
Magnus force field term B. This situation is the subject of the present 
article, in which we wish to show that the lifetime of the supercurrent 
can be calculated explicitely in terms of the phenomenological parameters of 
the material by means of a straightforward application of the field-theoretic 
relativistic formalism leading to the decay rate of the vacuum in scalar quan- 
tum electrodynamics. Indeed, the excitations of a relativistic quantum field 
encode in a rather natural way the particle-antiparticle quantum fluctuations 
of the vacuum and therefore the particle-antiparticle creation. This formal- 
ism, independently already proposed by Ping Ao |l3|], albeit for the simplest 
of the cases treated herewith, allows us to evaluate T for a number of vortex 
dynamics cases relevant to experimentally accessible situations. With refer- 
ence to the work of Ping Ao |l3|], it would appear that a main difference with 
the present work is the fact that the effects of the Magnus force, of dissipation 
and/or of pinning were not included in the resulting formula for the vortex 
pair production rate. For all the above-mentioned cases, the central technical 
ingredient of our formulation for the zero-temperature dynamics is an effective 
screened interaction and the resulting quadratic form of the effective-particle 
Lagrangian for the Feynman path-integral formulation of vacuum decay. 

Indeed, our article is organised as follows. In Section II we sketch some 
basic facts about the vortex dynamics and introduce our quantum field theory 
(QFT) formulation. This is expanded in Section III, where we treat in detail 
the situation in which dissipation and pinning are both absent, with the su- 
percurrent inducing a vortex-antivortex pair production in the presence of the 
homogeneous fields of the supercurrent and of the Magnus force. Two phe- 
nomenological parameters are introduced to describe the isolated vortex: its 
nucleation energy £q, suitably renormalised (e.g. to include screening effects 
by the external current or by the vortex plasma above the Kosterlitz-Thouless 
(KT) transition and pinning by impurities), and the vortex inertial mass m. 
We consider m and £q to be independent phenomenological parameters so 
as to account for the issue, debated in the recent vortex-dynamics literature 
PJll], where the inertial mass m may be taken as negligible (but see also Q). 
It is seen that the Magnus force raises a threshold for the magnitude of the 
supercurrent (or, the "electric" field) in which the pair production may occur. 
This was also qualitatively observed in Q|| ) . In Section IV we include the ef- 
fects of quantum dissipation through the formulation of Caldeira and Leggett 
p"2| , in which the drag force experienced by the moving vortex is mimicked by 
a linear coupling to the coordinates of an infinite set of harmonic oscillators. 
The main qualitative effect is seen to be the suppression of the threshold for 
pair production, as the Landau-level stability of the vortex in the field of the 
Magnus force is destroyed by any infinitesimal amount of dissipation. The 
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formula for the production rate T is worked out in detail for the dissipation- 
dominated case, that is in the limit 7 = m/£o^O. We deal mainly with the 
ohmic-case, to then show that sub- and super-ohmic cases can be treated in a 
similar way to yield a less pronounced pair-production phenomenon. In Sec- 
tion V we discuss the possibility of including the effects of pinning centres. A 
finite pinning barrier in the direction perpendicular to the current essentially 
results in a contribution to the activation energy £q, however pinning in the 
direction of the current has the consequence of suppressing the effect of the 
Magnus force in the dissipation-dominated dynamics. We also discuss the 
extreme case where pinning can be mimicked via an (unbounded) anisotropic 
harmonic force. For large enough separation, this corresponds to a confining 
linear potential in our relativistic formulation and this may mimick a situa- 
tion in which a uniform pinning pressure is exerted. The result is that, as 
soon as a component of the force in the direction perpendicular to the current 
is present, a threshold for vortex production appears. Section VI contains a 
more tentative discussion on how to implement the results obtained by the 
QFT approach in order to evaluate the vortex-induced voltage drop in the 
direction of the current across the sample, including also the effect of vortices 
which may be activated by thermal fluctuations. Finally, in the Appendix we 
show that some of the qualitative features observed within the QFT approach 
can also be seen to agree with a study of the classical equations of motion for 
a single vortex. 



II. VORTEX DYNAMICS: A QUANTUM FIELD THEORY APPROACH 

We begin this Section by briefly recalling the description of vortex dy- 
namics based on the time-dependent Landau-Ginzburg (TDLG) formulation 
for the superfluid system in the presence of vortex solutions and of a super- 
current. For more extensive discussions, also of related theoretical issues, the 



reader is referred for instance to the articles of the Seattle group || , |b|,|9|] and 
that of Lee and Fisher |L5| (see also [jlO|). For a superconductor, the TDLG 
effective action for the wavefunction ip(x, y, t) = y/pje^ 2 ^' 4 ) corresponds to 
the Lagrangian (in the units in which ft = 1 and the carriers' charge e* = 1): 

£ = -Ps0 - ^—p s (w - a) 2 + f( Ps ,Ps, vp s ) (2.1) 

Here uiq is the effective mass of the carriers, p s the 2-D superfluid density and 
/ contains the additional, phase independent, terms of the LG expansion in 
and its derivatives. In the following, we shall make the approximation of 
constant p s , thus dropping f(p s , p s , Vp s ), and treat the vortices as point-like 
objects. The external supercurrent is related to the external vector potential 
A via 

J = 7T = — (V0-A) (2.2) 
5 A niQ 
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Our situation consists in assuming that J is fixed by some external device 
"driving" a constant current which we take to be uniform throughout the 
sample. We are then interested in estimating the longitudinal potential drop 
and thus the ohmic resistance (in general current-dependent) due to the pro- 
duction and motion of the vortices. We implement the requirement of fixed 
current J by adding the term AJ to the Lagrangian and treating A as a 
Lagrange multiplier: 

C = -PsO - ^—Ps (V0 - A) 2 - A-J 

zmo 

Then 5C/5A = fixes the current: J = ^ (V0 — A); substituting back and 
disregarding a constant term -jr-J 2 , we get 



£ = - Ps 6 - V9-3 

A vortex is introduced as a configuration of the phase 6(r,t) having non-trivial 
topology and such that V# has formally a non-vanishing rotational (a scalar, 
in 2D): 

VxV6» = p v (r) 

In the presence of such singular vortex configurations with topological charges 
e n = 2im (where, in the following, n = ±1) and treating the vortices as 
point-like particles, we can rewrite the action in terms of a vortex density 
P«( r ) = En e nS(r - r n ) and vortex current density J v (r) = Y.n <5nr n J(r - r n ), 
with r n (t) the vortex position. The result is as follows 

S = J dtd 2 rC = J dtd 2 r (-p v V v - J v -A v ) 

J \ n n J 

where we have introduced the vortex fields A„ and V v through VxA„ = 
B = p s and E = —W v = xj. One can formally obtain the Magnus 
force term — J^-A^ from the — p s 9 term in the TDGL expansion (|2.l| ), in 
a shorthand way and disregarding possible subtle questions related to bound- 
ary conditions (for a more thorough approach see for instance [|| p^|JT5[| ). If 
#( r ) = En e n0 (r - r n (t)), then 9 = - J2 n e„r n -V6'(r - r n ). If we denote by 
A" 1 the inverse of the Laplacian operator, then 9 = e^ u A 1 d u J Vfl , where the 
vortex current is J Vfl = J2 n e nX n ^{^ — r n ) and we have written, from Eq. 
flU), d^{r - r n ) = -e^A-^Sir - r n ). Therefore 



J d 2 rp s 9 = J d 2 rA v -J t 



having formally put A Vfl = -e^A 1 d u p s . Also, it is straightforward to show 
that 



J d 2 rJ-V9 = - J d 2 rVV v -(xV9) = J d 2 rV v p v (21 

Adding a kinetic energy term J2n 5 mr n to the Lagrangian, Eq. (|2.6| ) shows 
that on the vortex acts a Lorentz-like force linked to a magnetic-like field B 



of intensity p s . We also recover the classical equation of motion (LI) for the 
vortex: the electric part of the field acting on the vortex is associated with the 
dual vector of the externally-applied supercurrent J. The picture is therefore 
quite reminiscent of the situation for scalar two-space-dimensional quantum 
electrodynamics (QED). 

In practice, the "electric field" E = xj will be really uniform only on 
some scale for which the current J is uniform. For scales less than that, the 
vortices themselves will alter the uniformity of the current. In the language 
of QED, a vortex at position r = ri near an antivortex at position r = r2 will 
also feel - beside the background electrostatic potential V v (ri) corresponding 
to the uniform electric field - the two-dimensional electrostatic attraction due 
to the antivortex 

AK (ri ) = ^ln^^i (2.9) 

where a is some scale which we can conveniently take to represent the vor- 
tex size. The above vortex-antivortex "electrostatic" interaction gives rise 



to an additional barrier, which can be taken into account [16 ,3] by an ad- 
ditive renormalization term in the activation energy. Following Minnhagen 
[16] (see also §), we make use of the fact that the barrier maximum occurs 
for |ri — ro|~-^ci = drS thus contributing a term to be included in the 

I x moil mo J ° 

renormalization of the activation energy of a vortex (e\ = 2ir) 

S OR = S -^ln(^) + ... (2.10) 



em V Ps 

Notice the appearance of the dielectrict constant e = e(T) which keeps into 
account the polarization due to the thermal effects below the KT transition. 
Above the KT transition, the "electrostatic" interaction is screened by the 
KT screening length Xkt and the renormalized activation energy will be, for 
J<^p s /m X KT 

£ 0R = £ -7r^ln{a/X KT ) + --- (2.11) 
m 

In the following we will lump these renormalization effects in the definition of 
the activation energy, which we call simply £q. 

The central issue of the present work is that the "electric" field term p v V v 
in Eq. ( |2.6| ) generates an instability in the system, which will have a tendency 
of producing vortex-antivortex pairs. We propose to describe pair creation 
through the formalism of relativistic QED. In this formalism, pairs are created 
out of the vacuum in the presence of a constant electromagnetic field, much as 
in our electromagnetic-analogue situation for the vortex dynamics. Further 
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terms will have to be added to the action ( |2.6| ), which describes only the 
interaction with the electromagnetic field, in order to account for the vortex 
motion. There will be, first of all, the energy Sq of the isolated vortex at rest in 
the absence of the external current. This we take as an effective renormalised 
activation energy Sqr to include the electrostatic screening energy AS as 
well as the average energy required to overcome possible pinning potential 
barriers. We consider, to keep the discussion as general as possible, also a 
kinetic energy term yr 2 , where m is the possible inertial mass of the vortex. 
In the relativistic formalism, which we employ, the total "free particle" energy 
for a vortex having momentum p would then be 



S = J£§ + -p 2 ^o + -Lp 2 + • • • = So + ^r 2 + • • • (2.12) 
y 7 2m 2 



Notice that in (2.12) the parameter - = So/m plays the role of the square 



of the speed of light c, in that the limit 7— >0 reproduces the non-relativistic 
(c^oo) expression for the kinetic energy. We imagine that the vortices will 
have a rather small kinetic energy as compared to £q. Indeed we will introduce 
the important effects of dissipation, by means of the Caldeira-Leggett formu- 
lation, and consider in particular the interesting case in which dissipation is 
more important than inertia ||. This formally corresponds to the limit m— >0 
or 7^0 at fixed rj. Note that in this case the drift velocity in the direction 
parallel to the electric field is, classically 

<*>■ " W&h < 2 ' 13 > 

Thus, we will implicitely look at the non-relativistic limit of a relativistic for- 
mulation. Notice that in the above formula and in the remainder of the paper 
we redefine the fields according to E^27rE and B^27rB, unless otherwise 
stated, so that our particles have formally charges e = dbl. 

The free particle relativistic expression of the energy, Eq. ( p. 12 ), corre- 
sponds to a space-time anisotropic scalar QFT with Lagrangian 

£ (<A) = <9o0*d O (/> - -V0*-V0 - 5 O 2 0> (2.14) 
7 

Notice that we have introduced a relativistic complex scalar field <f> describing 
both vortices and antivortices as its particle-antiparticle content. Thus, for 
a relativistic QFT in which particle pairs corresponding to vortices and an- 
tivortices are nucleated homogeneously from an external field A^=(V V ,A V ), 
the Lagrangian is given by 

C(<f>) = D (j)*D (j) - -D0*-D0 - Sl<f>*<j) (2.15) 
7 

in which we have introduced the covariant derivative = — iA^. There- 
fore, we see that the relativistic formulation captures quite naturally the quan- 
tum version of a plausible situation in which the activation energy and the 
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inertial mass are in general unrelated in the vortex motion through the su- 
percurrent. We are therefore in a position to evaluate the production rate T 
of the vort ex- ant i vortex pairs from a uniform "electromagnetic" field as a 
function of the field strengths E and B and of the ratio 7 = tji/Sq. 



III. AN EXACT FORMULA FOR THE VORTEX- ANTIVORTEX PRODUCTION 

RATE 

The calculation that follows is the two-dimensional scalar version of the 
well-known Schwinger calculation for the three-dimensional quantum electro- 
dynamic problem of vacuum decay by production of electron-positron pairs 



that are taken to infinity upon creation [17]. We evaluate the probability 
amplitude for the vacuum decay in time T 

Z = (0\e- lAT \0)=e- lTWo (3.1) 

where Wq = £(vac) — is taken to give the energy of the vacuum £(vac) 
and its decay rate T. With a suitable normalization factor, the probability 
amplitude is given by the functional integral over field configurations 

Z = M Jv<Pexp j-i J d 2 rdt(f>* (-D 2 , + ^D 2 - j (3.2) 

which is conveniently evaluated - formally - in the Euclidean metric 

Z = exp {-Trlii f--D 2 - D\ + £$ 
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exp(-*TW ) (3.3) 



by means of the known identity 

-D\ + £l ,. r fir 



Tr In B A T =lim/ ^Tr {^1+^ - e ~^} (3.4) 

in which we have written the square of the covariant anisotropic Euclidean 
Laplacian (with x 3 =it): D 2 E = D\ + ±(D 2 + D\). 

The evaluation of the trace is straightforward with the method of the 
Feynman path-integral; for the (d + l)-dimensional free particle the result is, 
for a system of size L d xT: 

Tr {e-^ T } = J dg (go|e- ( - 9|)T ko) 

= / ? (0)= 9 (.)= ?0 Pg( * )eXP {"/o Td ^ m ^^4 



iTL 



where we have employed a fictitious anisotropic mass term m^: mi = m% = % 
and ttt-3 = \. Hence Wq is real and, of course, there is no vacuum decay. In 
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the above manipulations q is taken to be a (d + l)-dimensional coordinate 
with anisotropic Euclidean metric; for the particle coupled to a gauge field 
case, we have the analogous path-integral construction 

Tr (e-(- D !H = [dq [ Vq{t)e~ & 

l ' J J q (o)=g(r)=q 

where, with the anisotropic fictitious mass terms, the Euclidean version of the 
relativistic Lagrangian for the vortex is 

The uniform-field situation corresponds to A^(q) = ^F^ v q u , with the (2+1)- 
dimensional field tensor given by 



/ B IE. 



1 [tv 



X 



-B iE y 
\-iE x -iE y 



after the analytic continuation to Euclidean time. The diagonalization of 
the operator (—D^)t is then achieved via Fourier transformation in the 
Euclidean-time variable: 



FT T 

/ dtC E = — V ttuv(0J n )qJUn)qv(-Un) 
JO 27T 



where 



4 w n 2 2 x 



%u(0On) = ~fB ^iE y 



J y 

2 

2 1 x 2 1 V 4 



The path-integral is now readily evaluated, to give 
-(-D|)r| =iTL 2 N{t) -Q j^-detol 

n=l ^ ' 



a-tWM n | ) V} 1 fl {i + i (~ + 5 



7r^n^ V 7 7 



where the normalization factor J\f(r) is adjusted so as to give the free-particle 
result, Eq. (|3.5|), in the limit of zero field. Using the known infinite-product 
formula n^=i(l + x 2 /TT 2 n 2 ) = sinh(x)/x, we see that we get 



r r {e-<-!>7 = s W-i-) 



r,/f-f 



V4W ^(ry/S^) 



At this point, the formula for the supercurrent decay rate becomes, from ( |3.3[ ), 
and (EUD 
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£ 2 B 2 



where a suitable vacuum subtraction has been inserted (so that no pair pro- 
duction takes place in the absence of an external field). This formula can be 
evaluated by summing over the residues at the poles of the or/sin (a;) func- 
tion, namely x=t^J~^- — ^ = nir, except for the pole at x = 0. To give 
meaning to the poles one attributes a small negative imaginary part to x, so 
that Irn-^r- = fm = (— ) n ir5(x — irn). The end result is, for the pair 

production rate per unit area 

The main result in this exact expression, valid for any ratio 7, is the presence of 
a threshold for the decay of the supercurrent, E > E C =B / This conclusion 
could have been obtained from the classical equation of motion, Eq. ( |1.1|) . 
In fact, from ( |1.1| ) the classical trajectory, without dissipation i.e. r] = 0, 
is a cycloid describing a motion on average orthogonal to E, that is with 
the vortices dragged along with the current J (see Appendix). The average 
displacement in the direction parallel to E turns out to be (Ax) = Em/B 2 , 
thus we see that a vortex can nucleate if the energy gained from the electric 
field can at least equal the nucleation energy: E(Ax)>£q, which gives the 
above threshold condition. However, even though vortices can nucleate above 
threshold, the fact that there is no net transport in the direction orthogonal 
to the current in the absence of dissipation implies that no voltage drop will 
result in this case, hence no resistence. 

In the following we will discuss the interesting and realistic case in which 
there is always friction, induced by the quantum dissipation (for a discussion 



of a microscopic theory see [14]). In this case particles will drift, classically 
and after a transient during which a few cycloidal spirals are noted, along 
a straight line with drift velocity in the direction of E given by Eq. (2.13). 



Therefore we anticipate the absence of a threshold for this case, and since a 
net transport in the direction of E is activated we conclude that a voltage 
drop will ensue. 

As a final comment on the Eq. ( 3.14j ) we note that, in the absence of the 



Magnus field B and at the n = 1 order - also remembering that below the KT 



transition £q depends on J as in Eq. ( 2.10 ) - the argument of the exponential 



has the same leading dependence on the current J as reported in the work of 
Ping Ao |13|. 
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IV. VORTEX NUCLEATION IN THE PRESENCE OF QUANTUM 

DISSIPATION 



Much work in the physics of vortex dynamics assumes the presence of 
an amount of dissipation. At a quantum level, this can be treated with the 



formalism introduced by Caldeira and Leggett [12]. This assumes the moving 
quantum particle to interact with a bath of quantum harmonic oscillators 
having an arbitrary distribution of masses m k and frequencies uj k , the only 
physical constraints being that the coupling between particle coordinate q 
and harmonic coordinates x k is linear and the classical equation of motion 
reproduces the form ( |1 . 1[ ) . For the problem at hand, therefore, the Lagrangian 
describing vortex motion in the presence of a supercurrent and dissipation is, 
in Euclidean time: 

11 f 1 c 2 

where, as is well known, the last term in the oscillators Lagrangian is added 
to avoid the lowering of the minimum in the particle potential by means of 



the coupling to the external bath, and where the c k are constrained by \12\ 



2 m k u k 

r\ being the phenomenological friction coefficient. The above constraint cor- 
responds to the so-called ohmic case for the dissipation; more generally, the 
spectral function can have a frequency-dependent form {lo c is a cutoff fre- 
quency) J(u) = rjuj s exp(— uj/oj c ), where s > 1 corresponds to the super-ohmic 
and s < 1 to the sub-ohmic cases. Although we deal mainly with the ohmic, 
s = 1 case in what follows we will also extend our results in a rather simple 
fashion to the non-ohmic cases in order to show how vortex-pair production 
is optimised in the ohmic case. After Fourier transformation, the x^-modes 
can be integrated out, for example by means of their equation of motion in 
terms of the action = J Q T dt£g 

dS D t f \ 
Q^fy = ^ (™^ 2 + w2)*fc(-") + c fe q(-^)J = 

leading to the effective action 

S E = \qfiM%uMqu(-<^n) 7, 57-lh 2T<7a(Wn)<?a(-W n .)} 

where the oscillator part contains a summation over the space- like components 
a=l, 2 only of the particle coordinate. The sum over uj k can be carried out 
by resorting to the constraint fl4.2j); in fact 



c l f + °° doo ^ c| _ 7} 

1^ 0™,,.,2A,2 T ,.,2N - /„ ,.,2 , ,.,2 2^ 77— ~2 d ^ ~ U k ) 



2m k u%(u^ + u}%) Jo u>n + uj 2 ^ 2m k u k " 2\io n \ 
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We therefore attain the form 



2vr ?' 



^(w n )^(o; n )^(-w r , 



with Q D modified only in its diagonal space-like elements by the inclusion of 
dissipation; for u n > 



/ j, ,2 , 2, , 



V 



2 L 
2 

" ~^^E X 



2 - 
4 W n + 2 W ™ 
2 '-^2/ 



2 

,.,2 



Notice that this implies that in Euclidean time dissipation processes can be 



represented by means of an effective action. At this point, the formula ( 3.13 ) 
for the decay rate of the supercurrent, in the presence of both the uniform 
external fields and of dissipation, makes use of the matrix determinant 



det Q D 



T 



7 2 7 7 

— id + -U) 

4 2 



T 



4 2 



16 



E 2 + ^B 2 



which, adopting a factorization similar to the one employed without dissipa- 
tion, leads to 



/oo A I oo 

— iM'(T)e- £ o T { n 
T n=l 



E 2 



+ 



uj„B 



r]u n /2 + 7^/4 4{ V uj n /2 + 7<4/4) s 
with the normalization factor JV'(t) given by 



AA(r)I] 



n=l 




*7 7 2 



7 



\Aitt 



3/2 



n 



71=1 



27] 



Notice that we have explicitely rewritten the vacuum subtraction term in such 
a way that no pair production takes place in the absence of an external field, 
E = B = 0. In this case a closed-form expression for V does not appear 
possible; nevertheless the leading contribution to the r-integral can be found 
from the residue of the main, n = 1 pole of its argument. No relevant pole 
arises from M' (r), hence the main contribution is from w = u\ = 2it/ti which 
is the real positive zero of 



1 



E 2 



+ 



uj 2 B 2 



r]u /2 + 7w 2 /4 A(t]lo /2 + 7w 2 /4) 2 







For given B, rj and 7, this cubic equation has a real zero for any infinitesimal 
electric field E, as indeed for E— >0 we get 



2rj 



7] 2 + B 2 



E + 



The physical meaning of this is that while in the absence of dissipation, with 

IE 2 B 2 
7 T 5 "' 



real root appears only for E > B/y^y, now any in- 
finitesimal E, therefore any infinitesimal supercurrent, gives rise to vortex 
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pair production and decay. The main consequence of (ohmic) dissipation is 
therefore the elimination of the threshold for a non- vanishing T, as was also 
noted by , e.g., Ao and Thouless M. Noticing that uj\ is independent of 7, we 
work out the final formula in the "dissipaton-dominated" case, that is in the 
limit 7— >0 also considered by Stephen || , although a less instructive analytic 
result could be obtained for arbitrary 7 to leading order in E 2 . The formula 
for r is first of all rewritten, more conveniently, as 



2L 2 



Re 



dr. / 1 \ 



"'7 v a 

T \ 4-7TT 



\ 3/2 00 ( 

n 1 

7 n=l I 



+ 4 



B 2 + r] 2 + j(r]U n - £ 2 ) 



7 2 <4 



n 1 



8t]E 2 



-1 



n=l 



7 2 o;3 + 4 Wn (£2 + r/ 2 + 7(770^ - E 2 )) 
which for 7^0 simplifies to 



lim 



2L 2 



Re 



2r/£ 2 



UJ n (B 2 + T] 2 ) 



with the following asymptotic form of the normalization factor 



M"{r) 



7 



/ 1 \ 3/2 00 / 
n=l \ 

Attt J 7 v 



B 2 + 7? 2 



in which the limit form xj sinh(x)— >2xe~ x has been used. The leading pole 
has the precise form (4.12), but for arbitrary E, and evaluating the residue 
at this pole the formula for T reads, with £ 2 R = £ q + -\/B 2 + r/ 2 



lim 



2L 2 



47TT1 



3/2 



£ 2 + r? 2 IJ 1 



ra=2 



-1 



The last infinite product is clearly divergent unless a frequency cutoff is in- 
troduced, id n <u> c where uj c is a large frequency above which the effects of dis- 
sipation can be neglected [jl2]. This implies the existence of a cutoff integer 
n* = [u c /ui] and the last factor will thus be a cutoff-dependent dimensionless 
number C(n*) = C* . For small E 2 we see that n* = is large and in 

this case we can estimate the asymptotic behaviour of C* by 



n=2 



n 



■ exp — n 



cixln f 1 



2/n* 



1 



e * 
4 



67/ 



r? 2 + S 2 2^ 2 



If, on the other hand, E 2 is very large and formally n* shrinks to zero, it 
would mean that our situation is such that the Caldeira-Leggett dissipative 
regime breaks down. We assume here to be always in a range of E compatible 
with the Caldeira-Leggett regime, thus n* can be either large (corresponding 
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to C* given by ( 4.17| )) or finite (corresponding to C* of order unity). The final 



asymptotic form for the pair production rate when 7— >0 will be, therefore, 
substituting t\ = 2tt/loi from Eq. ( |4.12|) in Eq. (|4.16|) , 



where we have defined the "renormalized" activation energy 

£ 2 0R = £ 2 + ^B 2 + r ] 2 (4.19) 

This has a physical interpretation connected with the energy B/2m of the 
lowest Landau level which becomes infinite in the limit of zero inertial mass 
and contributes a diverging zero-point energy to the total activation energy. 
Indeed, for 77 = 0, £qr— £q + B 2 /2m. If an infinite barrier for placing a vortex 
in the lowest Landau level does really develop, then vortex pair production 
will cease and there will be no supercurrent decay. We therefore generically 
assume in the following that £qr is ultimately just a phenomenological finite- 
value parameter of the theory. 

At this point we can generalise the results obtained above for the ohmic, 
s = 1, case to the non-ohmic, s / 1, situations (for a recent discussion on 
the relevance of non-ohmic cases, see [pi ). Indeed, Eq. ( 4.12| ) holds good 



if J(co) = r]u) s provided we replace r] with a frequency-dependent friction 
coefficient r/(uj) = r]uj p , with p = s — 1. Then we have the following equation 
for the main pole tu\ = 2h/t\ 



r\ A io^ + B A 



which we can qualitatively handle in the limit E— »0. First, we observe from 
Eq. (4.20) that the role of B is enhanced in the super-ohmic case, while it 



is suppressed if dissipation is sub-ohmic, in agreement with ||]. Then, both 
in the super-ohmic p > and in the sub-ohmic p < cases, we get that 
o;i~£' 2 /( 1 ~l p l) so that in the relevant formula for T, e.g. Eq. ( [4.18] ), we must 
substitute E with £ ,1 /( 1_ IpD. For small E this leads to a suppression of the 
vortex production rate as soon as p ^ 0, compared to the ohmic case which 
therefore represents the situation with the maximum production rate. 

To conclude this Section, we point out that the argument of the exponen- 
tial appearing in Eq. ( (4.18|) has a direct physical interpretation. Apart from 
constants, it is the Action (=EnergyxTime) required for the nucleation of a 
vortex, since Energy=£oR and Time=^Ar /vn = ^^'^e • Indeed, In = is 
the distance involved for the nucleation in a field E whereas vn = is the 
drift velocity (if we indicate rj{B) = {rf + B 2 )/r]). 
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VORTEX NUCLEATION IN THE PRESENCE OF A PINNING POTENTIAL 



Most superconducting materials in which vortex dynamics plays a role 
contain a finite density of defects of some sort. Many discussions in the liter- 
ature deal with the tunneling effect of vortices trapped by a potential barrier 
on the resistence of a superconductor. A finite height and width of the pinning 
barrier in the direction of the "electric field" E, that is perpendicular to the 
current J, ultimately give rise to a finite additional term which renormalises 
(?o_r- Pinning barriers in the direction of the current have instead a different 
physical effect, in particular in the dissipation-dominated case. In our own 
relativistic formulation we can account for an additional effect due to pinning 
if we consider the "extreme" case of the unbounded harmonic-well potential 

U( q ) = ~k x ql+~k y q 2 y (5.1) 

which we add to our quantum field Lagrangian through the term U(r)4>*(fi. 
Strictly-speaking, in our relativistic treatment this amounts to a confinement 
force that, while linear in q at small distances, at large distances remains con- 
stant and thus simulates the presence of a uniform "confinement pressure" on 
the vortices. Thus, in this "extreme" case a threshold for the pair production 
may reappear. In what follows, we first treat the case of the unbounded po- 
tential in the absence of dissipation, and then discuss the more realistic case 
of a finite potential barrier. Technically, with the form ( |5.1[ ) we can formally 
repeat the calculation of the previous Sections, with a quadratic form in the 
path integral defined by the matrix 



/ lu 2 n + %u n + \k x ^iE x 

-Tf£ l^ + R^ + lky ^iE y | (5.2) 

\ 2 ^ 2 y i 



Having the limit 7^0 always in mind, we begin the discussion of this ex- 
treme case of pinning by considering the absence of dissipation, r] = 0. Then, 
with a suitable factorization for det f2 we can write, for the pair production 
rate and in the limit 7^0 

T f°° dr , & f -S- / k x k v - 2{k x El + k v E*)Y l ) 

lim — =ftejf ^ iA r(r)e^| n( 1 + ) " const - | ( 5 - 3 ) 



The normalization factor has again a singular limit form 



n=l 

/ 1 \ 3/2 T B 

in which we see that the very same quantum zero-point energy fluctuations 
in the Landau levels give rise to a divergent renormalization in the rest mass 
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for 7— >0: £q^£q R = £q + B/7. The remainder of the integral in r can be 
evaluated as in Section III in this limit, and we end up with the exact formula 



lim F = "2^ (^ 2 + ^ " ^^Z 2 ) } 

. £ (_l)n+l n -l/2 exp J ^Ofi_ 



n=l 



\{k x El + k y El-k x k y /2) 



-.n 



in which we see that Y is zero unless the condition 2E X /k x + 2i? 2 / > 1 is 
satisfied, that is no production of vortex-antivortex pairs takes place inside 



an ellipse in (E x ,E y ) having half-axes \Jk x j2 and Jk y /2. The physical in- 
terpretation of this phenomenon is that the threshold in the minimum value 
of the electric field for a non-vanishing pair production is now restored by 
confinement. Indeed, while particles get trapped in stable Landau levels for a 
vanishing inertial mass, under the action of a pinning potential the trapping 
role of the magnetic field is restricted to a renormalization of the rest mass. 
Transitions out of these levels can occur for strong enough electric fields since 
the situation becomes now a quasi-static equilibrium between the confine- 
ment and the electrostatic forces, the dynamics being separated away with 
the magnetic field term. Notice, however, that for E y = and k x = (current 
and confinement both in the y-direction) there is a non-vanishing rate T for 
any E x ^0, hence there is no threshold. Nevertheless, no net transport takes 
place in the direction of the electric field and the supercurrent does not decay 
if dissipation is absent (see Appendix). 

We now consider the case in which ohmic dissipation is present in a more 
realistic setting, by specializing to E y = k x = 0. In fact, one can take into 
account a finite-width pinning barrier in the x-direction by a renormalization 
of £q, but a new qualitative feature appears if pinning is in the direction of 
the current, J = — xE. We adopt the following factorization for the matrix 
determinant: 

det dP = ^ f V) 2 { 1 + V + + + ^ ~ ^ 




7 2 u> 2 



(4k y - 8E 2 ) t]oj- 1 - 8k y E 2 uj 



2, ,-2 



\ 7 2 w 2 + V + 4B 2 + 7(4770; + 2k y - 4E 2 ) 
since for 7— >0 this becomes 

limdetO- = ^fVV(l + 4^!} 

v + — y — 




B 2 + 7] 2 

Here we see that the prefactor is of the type already seen for dissipation, 



hence it gives rise to a nucleation energy renormalization of the type (4.19). 
Denoting 
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B 2 = B 2 + rj 2 

E 2 =E 2 ~ hy 

K 2 = h y E 2 (5.8) 

we see that the dominant pole in the integral formula for V is the solution of 
B 2 vo 2 — 2r]E 2 u>i — \K 2 = 0. The only relevant principal pole is then for 

U! = y.= (^E 2 + ^ri 2 E* + 4B 2 K 2 ^ /B 2 (5.9) 

for which the residue is readily calculated and yields the main contribution 
to the vacuum decay rate 

r / i \ 3 / 2 / 

lim—^2TTe- £ oR Tl i^—J C*^B 2 + r] 2 (5.10) 

where C* is expressed by an infinite product, cutoff at n* = u c /uj\, similarly 
to what seen in Section IV. The formula for the pole, Eq. ( |5.9| ), reveals its 



physical contents in the limit in which the E 2 <^k y . Expanding the square-root 
to lowest order, we always obtain a pole for small E 

E 2 

wi = 2 — (5.11) 
V 



which is of the same form as Eq. ( 4.12] ) for the purely dissipative case, except 



that - due to pinning - the renormalization effect induced by the Magnus field 
on the friction coefficient has disappeared. This result is in agreement with the 
analysis based on the classical equations presented in the Appendix, which, 
furthermore, shows that Eq. (5.11) holds (for small E) also for a finite-height 



and finite-width pinning barrier. In the general case in which the pinning cen- 
tres will be distributed with some density, this will correspond to an effective 
friction coefficient that will present a dependence on B intermediate between 
the two cases: r\ < r) e ff < {rj 2 + B 2 )/rj. 



VI. AN APPROACH TO THERMAL EFFECTS AND CONCLUSIONS 

To complete the picture one would need to evaluate the vortex nucleation 
rate under the effect of both the external current J (resulting in the "electric" 
field E) and the temperature T in order to estimate the voltage drop. This 
appears to be a rather difficult task, since it would involve a theoretical de- 
scription of how the thermal distribution of vortices is obtained from a state 
with no vortices, a process of non-equilibrium thermodynamics. 

Here we will take a less ambitious view and try to obtain a phenomeno- 
logical description of the process of vortex pair production in the presence 
of thermal fluctuations. By a detailed balance argument, to be developed 
below, one finds that the density of free vortices pf is proportional to the 
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square root of the production rate, pf~y/T. Considering the relevant case of 
the dissipation-dominated dynamics described in Section III, we tentatively 
propose to simulate the effect of temperature on the vortex production by 
introducing thermal currents J(T) (induced for instance by phonons, or oth- 
erwise) and using the result of Eq. ([OJ|), with E 2 = (2irJ(T)) 2 , namely 
r~exp —SojiVeff /47rJ 2 (T). The requirement of obtaining a Boltzmann dis- 
tribution for pf fixes J 2 (T) = ^7] e ff£onT (we take kg = 1 for Boltzmann's 
constant in what follows). In the general case, where we have both tempera- 
ture T and external current J, we will use Eq. ( |4.1g| ) with E 2 = {2tt Jtot) 2 by 
adding incoherently the two contributions: J 2 ot = J{T) 2 + J 2 . We can thus 
rewrite, in the general case 

' C*^jJ tot exp -2£ 0R /(T + AT(J)) 



L X Ly 



where L XtV are the linear sizes of the sample; also 
AT(J) = ~ UJ)- 1 

and we have indicated that AT (J) is proportional to the inverse of the nucle- 
ation time for a vortex (see the discussion at the end of Section IV). Let us 
recall that below the KT transition £qr depends on J (see Eq. ( |2.10|) ) and in 
this dependence we include just the external coherent current as < Jtot >= J- 
Therefore, vortices are not produced for zero current and temperatures below 
the KT transition. 

We can then compute the free vortex number density pf by a detailed 
balance argument: the production rate must equal the annihilation rate. Ne- 
glecting the rate at which the vortices disappear from the sample due to the 



drift described by Eq. ( 2.13 ), thus assuming macroscopic dimensions for the 
sample, the annihilation rate per unit area will be due to two mechanisms: 
1) A possible vortex-antivortex annihilation, and this rate will be proportional 
to the product of an annihilation length a (the 2-D analogue of the familiar 
3-D cross section), times the density of the free vortices pf, times the incident 
flux pf-2 < v >, where the average drift velocity will be given, according to 
Q2.13[ ), by (v) = 2tt J tot /rjeff- Hence: 



^^2p 2 (v)a 



2) A possible annihilation of vortices with pinned antivortices, with a rate 

ip f p p (v)a 



L X Ly 



where p p is the density of pinned vortices. In turn, at equilibrium we will 
have pjTp = p p T u , where T p and T u are the rates for pinning and unpinning, 
respectively. In conclusion, the detailed balance will be 
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L X Ly 



2 ^Jtot rn -p /T -, n 
Veff 



which gives 



Pf 



2vra(2 + T p /T u ) 



S { 



exp 



OR 



T + AT(J) 



Thus, the vortex current J v in the x-direction orthogonal to the supercurrent 
density J will be J v = 2irpfJ/rj e ff, yielding the potential drop Ay = L X J V 
in the y-direction parallel to J. Thus: 



AV = 2ttL^J 



C*rj 



Veff<r(2 + r p /r u ) 



exp 



T + AT(J) 



For small external currents such that AT(J)<^.T one gets, above the KT 
transition, in terms of the current / = L V J 



AV~R(T)I + Q(T)/ 3 + 



where 



R(T) = 2vr 



Q(T) 



R{T) 



rj eff a(2 + T p /T u ) 
8vr 



-£or/t 



Veff 



T 2 



It is to be stressed at this point that the above considerations have led us to a 
rather predictable result: there is no resistance in the sample if dissipation is 
absent. As we have previously pointed out, this qualitative conclusion could 
have been obtained from the classical equations of motion, Eq.s (|1,1| ) and 
( 2.13[ ) (see Appendix). We recall that, as discussed at the end of Section V 
and in the Appendix, in the general case of a distribution of pinning barriers 
one will have rj < rj e ff < (r/ 2 + 4-7r 2 p 2 )/r/ (remembering that B = 2irp s ). 

We finally point out that, for temperatures below the KT transition, £oi? 
depends on J, as in Eq. ( 2.10 ). Thus, for low J such that AT(J)<CT, 
we recover the standard result that AV~I 1+a , where a^2 for T^T KT HI 



Our formula includes dynamical effects due to current-induced nucleation. In 
particular, by expanding for low J, we can repeat the above calculations to 
arrive at the following dependence on J 



AV~ci/ 1+a + c 2 I 



r3+a 



hi 



mo J 
p s a 



+ ... 



where c\ and C2 have expressions similar to those reported in ( |6.9[ ). 

To conclude, we have presented a plausible physical scenario for the ap- 
plicability of our exact formalism giving the decay rate of a superconducting 
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current and thus the resistence, or voltage drop, observable in a real sam- 
ple. Thermal effects have been accounted for phenomenologically through 
thermally-induced currents, and pinning forces by means of a harmonic po- 
tential. We find that in some extremal cases there is a threshold for the 
production of vortex-antivortex pairs, in that these are trapped in Landau 
or harmonic-oscillator levels, but that the threshold is eliminated, if extremal 
pinning is absent, by the addition of any infinitesimal amount of dissipation. 
In the case of finite pinning, or even infinite-pinning barrier in the direction of 
the current, the threshold disappears. Pinning barriers in the direction of the 
current suppress the effect of the Magnus force (confined to a renormalization 
of the vortex nucleation energy). The limit of a vanishing inertial mass has 
been shown to lead to considerable simplification in the algebra, as well as 
to new physical effects like the divergent renormalization of the vortex nucle- 
ation energy, resulting in the inhibition of vortex production for all the cases 
considered when the inertial mass is strictly zero. Further work is needed 
to account for the thermal effects in the nucleation of vortices within a fully 
microscopic theoretical approach. 
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APPENDIX 



Here we summarize some relevant results from the solution of the classical 
equations of motion for the vortex dynamics, Eq. ( |1.1| ) (for charges with 
e 2 = 1). We begin by recalling that in the presence of the fields E = Ex and 
B = Bz the solution corresponding to the initial conditions x(0) = y(0) = 
and x(0) = y(0) = is of the form 

mE 

= Jjpt 1 ~ COSUJ c t) 

/ \ mE E , , _ . 

y(t) = -—smuj c t + -t (A.l) 

which represents a cycloid with frequency u c = eB/m. This shows that the 
average displacement along the electric field's x-direction is fixed at {x(t)) = 
xq = mE/eB 2 so that there is no transport in the direction orthogonal to 
the current, hence no voltage drop. Also, we require the field to be such that 
Exq>£q to extract the particle from its classical trajectory, in agreement with 
our conclusion on the existence of a threshold in Section III. If on the contrary 
friction is switched on, the particle will drift after a transient along a straight 
line with a drift velocity in the x-direction given precisely by Eq. (2.13); 
in this case transport is activated and a voltage drop ensues. Indeed, with 
the same initial conditions as above, the motion in the presence of friction is 
represented by 

x(t) = J^-^ ( e -^/ m )* cos^t - 2a) - cos2a| + ff^ t 

rj z + B z y > rj + B z 

V(t) = Prl^ f e- (r?/m)t sin(w c t - 2a) + sin 2a| + ^—^ t (A.2) 

rj z + B z i. > rj z -\- B* 

where tana = rj/B, showing that (after a transient damped-cycloidal motion 
that can be neglected for t>m/rj) we have an average x-component of the 
velocity {x(t)) = eE/rj(B) with rj(B) = {rj 2 + B 2 )/r] as in Eq. ( |2.13|) . Also, 
there is no threshold value of E for the vortex-antivortex nucleation. 

Next we consider what might be, classically, the main effect of pinning on 
the vortex motion. We introduce a pinning potential approximately described 
by U(x,y) = \k x (x-x p ) 2 + \ky(y-y p ) 2 for \x-x p \ < i x and \y—y p \ < i y and 
zero otherwise, where (x p , y p ) are the coordinates of the pinning centre. We 
imagine that there will be many pinning centres in the sample, but discuss 
the local dynamics of a vortex around one of the impurity sites. Our particles, 
i.e. the vortices, will feel a potential barrier around (x p ,y p ). Like in Section 
V, we assume that the potential barrier in the direction of the electric field 
E = Ex will have the main effect of increasing the activation energy Sqr, thus 
renormalising it by way of the addition of a further term. We can thus treat 
this part of the problem phenomenologically. The barrier in the y-direction 
has a different implication for the dynamics. Consider first the frictionless 
case, T) = 0. The general solution of the classical equation, including the force 
due to U (y), is 
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u c 1 k y eE 2 /c„ 

x(t)-x p = Acoslut + - / t +-^yoT + x 

lo 2 mfey + B z eB 

EB 

V(t) -V P = -Asinur + — - b2 t + yo (A.3) 

where r = t — to and = ^Juj 2 + k y /m. We see that, after averaging over the 
oscillations and with general initial conditions, (y(t) — y p ) = m! ^+ B -2 t + yo 

and (x(t)-Xp) = xq + ^yo T + \ m^+B^ ^ ' so that under the effect of pinning 
in the y-direction (x(t) — x p ) can reach, after a suitable time, a large enough 
value and allow the vortex to gain nucleation energy from the electric field. 
Depending on the values of the parameters this can occur within the pinning 
potential width t y , in which case the pinning in the y-direction removes the 
threshold for vortex nucleation which occurred without friction and pinning. 
However, the classical equations indicate that even in the present case there 
is no net transport in the x-direction if there is no friction. In fact, the 
previous solution holds up to \ (y(t) — y p )\ r ^£ y , beyond which either there is no 
longer any pinning, and thus the motion will fall back on the familiar cycloid 
(resulting in no transport), or else the vortex will feel the effect of another 
pinning centre at (x' p ,y' p ). In the latter case we will have, at the initial time 
to of the new piece of the trajectory, (y(io) — y' p ) = —t' y . Therefore, at the 
later time t\ for which (y(t\) — y p ) = l' y (that is as the vortex moves through 
the y-range of the new pinning centre from y = y p — l' y to y = y' p + £' p ), we 
find that (x(t\) — x') = (x(to) — x' p ). The conclusion is that there has been no 
net transport in the x-direction. The situation is completely different if there 
is friction, rj / 0, as can be easily seen by solving the classical equations by 
formally taking mix = and my = for long times and adding the pinning 
force due to U(y). The solution is that, for long times (t > m/r/), {x(t)) = 
eE/r/ whilst {y(t) — y') = BE/k y m Thus, there is always the possibility 
of nucleation and transport. Comparing with the case without pinning, we 
see that the drift velocity in the x-direction is now determined by the actual 
friction coefficient rj < r]{B) = {rj 2 + B 2 )/i], whilst in y one has simply a 
displacement (thus, for low E at least, y(t) will remain in the range where 
U (y) 7^ 0). In the general case, where there can be pinning centers distributed 
with some density, the effective friction coefficient rj e ff will be somehow in 
between the extreme cases, t] < i] e ff < (rj 2 + B 2 )/rj. 
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